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Formation of 77-mesic nucleus, a bound state of an 77 meson in a nu- 
cleus, is reviewed in this paper. Three different theoretical approaches are 
used to calculate the binding energies and widths of such nuclei. The effect 
of ry-mesic nucleus in pion double-charge-exchange reaction is discussed. 
Experimental efforts by different groups to detect the nucleus are also dis- 
cussed. The ramifications of the theoretical and experimental studies of 
the bound state of in a nucleus are pointed out. 

PACS numbers: 24.60.D, 13.75.G, 21.10.D 

1. Introduction 

In the past, bound systems of strongly interacting particles in a nucleus 
or atom, such as hyperons and mesons, have provided valuable information 
about various aspects of hadron-nucleon interaction in a many-body envi- 
ronment. For example, experimental and theoretical studies of hypernuclei 
led to impressive advances in our knowledge of the AA'^ and TiN interactions. 

Until recently it was believed that the 77 meson plays almost no role in 
nuclear physics because the r]NN coupling constant is very small compared 
to ttNN and vrA^A coupling constants. This is in sharp contrast to vr- 
nucleus interaction which has been studied extensively. The situation has, 
however, changed recently in medium- and high-energy nuclear reactions 
where significant amount of pion induced r] production has been observed 
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at pion energies near 500 MeV [Ij . In view of this, there has been a surge of 
interest in studying the ry meson within the framework of nuclear physics. 
Some of the studies have led to interesting surprises, such as the prediction 
of the existence of a bound state of rj meson in a nucleus, termed r/-mesic 
nucleus, and the phenomenon of "mesonic compound nucleus" in high en- 
ergy pion double-charge-exchange reactions. In this paper, we will describe 
the formation of the ry-mesic nucleus and point out the relevant physics that 
we can hope to learn by using the rj meson as a nuclear probe. 



2. Formation of Eta-Mesic Nucleus 

The existence of r/-mesic nucleus was first predicted by us in 1986 [2j. It 
is a consequence of the attractive interaction between the rj meson and all 
the nucleons in the nucleus. The attractive nature of the interaction follows 
from the work of Bhalerao and Liu [3j who found, from a detailed coupled- 
channel analysis of vrA^ — > vrA^, vrA^ ttttN, and vrA^ rjN reactions, that 
near-threshold rjN interaction is attractive. 

The binding energy and width of the r/-mesic nucleus are calcu- 
lated by solving the momentum-space relativistic three-dimensional integral 
equation 

^ V'(k') + / dk < k' I y I k > VXk) = EtPik') , (1) 

using the inverse-iteration method of Kwon and Tabakin [3]. Here < k' | 
y I k > are momentum-space matrix elements of the ??-nucleus optical 
potential V, with k and k' denoting, respectively, the initial and final r/- 
nucleus relative momenta. The fi is the reduced mass of the ry-nucleus 
system and E = e^^ -|-ir^/2 = k? is the complex eigenenergy. For bound 
states, both e^^ and F^ are negative. Three different theoretical approaches 
to V are used to calculate and they are described below. 



2.1. Covariant r]-Nucleus Optical Potential 

The first-order microscopic ?7-nucleus optical potential, after using the 
covariant reduction scheme of Celenza et al. [5], has the form[2l[6] 

< k' I y I k > = ^ 1 dQ < k', -(k' + Q) I t(^)^jV->r,iV I 
j 

X I k, -(k + Q) > <A*(-k' - Q)<P,i-k - Q) , (2) 

where the off-shell t]N interaction tr^N^r^N is weighted by the product of the 
nuclear wave functions 4'*j(l}j corresponding to having the nucleon j at the 
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Table 1. Binding energies and half- widths (both in MeV) of ?7-mesic nuclei given 
by the full off-shell calculation. No bound state solutions were found for A < 10. 



Nucleus 


Orbital {n£) 




12q 


Is 


-(1.19-F3.67i) 


160 


Is 


-(3.45-F5.38i) 


26Mg 


Is 


-(6.39-F6.60i) 


40Ca 


Is 


-(8.91 -F6.80i) 


90Zr 


Is 


-(14.80 -F8.87i) 




Ip 


-(4.75-F6.70i) 


208p]^ 


Is 


-(18.46 + lO.lli) 




2s 


-(2.37-F5.82i) 




Ip 


-(12.28 -F9.28i) 




Id 


-(3.99-F6.90i) 



momenta — (k-|-Q) and — (k'-|-Q) before and after the collision, respectively. 
The ry-nucleus interaction V is related to the elementary t^A^ process by the 
kinematical transformations of Liu and Shakin [7] . The ijN invariant mass 
in the cm. frame of the rj and the nucleon j is given by 



= [{w - EcAQ)}' - Q'] 



+ Mn- 



< {Mr, + MNf 



( M^ + Ma 

2Mc,j \Mr, + MN 



(3) 



where Q, Ecj and Mcj are, respectively, the momentum, total energy, 
and mass of the core nucleus arising from removing a nucleon j of momen- 
tum — (k -|- Q) and binding energy | €j \ from the target nucleus having 
the momentum — k. Calculation of V involves full off-shell kinematics and 
integration over the Fermi motion variable Q. It requires knowledge of the 
basic triN^TjN at subthreshold energies [3]. All kinematic quantities are cal- 
culated using well-established Lorentz transformations. For near threshold 
r]N interaction, only one resonance of the A^* isobar has to be considered 
for each partial wave. They are ^n, Pn, and -D13 resonances. The rjN 
interaction parameters are taken from ref . [3j . Details of the calculation can 
be found in refs. [21 [6]. 

The results of the calculation are presented in table [TJ The nuclear 
wave functions in eq.([2]) are derived from experimental form factors with the 
proton finite size corrected for. From the table 1, it can be seen that r] can be 
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bound in nuclei with mass number A > 10. The ry-nucleus interaction is not 
strong enough to have a bound state in hghter nuclei. This is because there 
is a reduction in the strength of the ry-nucleus interaction at subthreshold 
energies. The number of nuclear orbitals in which rj can be bound increases 
with increasing mass number. 

It should be mentioned that the calculations do not include the effects 
of Pauli blocking. However, estimates of the blocking using local-density 
approximation show about 5% reduction in the widths [2j. The calculated 
widths of r/-mesic nuclei vary from 10 to 13 MeV in lighter nuclei and are 
between 15 and 20 MeV for the heavier ones. About 95% of the width is 
due to the decay of the ?7-mesic nucleus with the emission of one pion and 
one nucleon. The remaining 5% is due to the emission of two pions [8]. 
As absorptions of rjN N and r]NN NN are strongly suppressed 
kinematically, their contributions to the width are insignificant. Widths of 
the order of keV may be possible if ijNNN — > NNN can take place. In 
this latter case, the r] can then share its energy and momentum with three 
nucleons. 

The binding energy is very sensitive to the rjN interaction coupling con- 
stant g^j^j^, ■ The value of g^^j^* determined in ref. [3j is 0.77. This leads 
to a value of 0.28 fm for the real part of the s-wave rjN scattering length 
0"riN, corresponding to an attractive interaction. As an example, it has been 
shown by Haider and Liu [9j that the bound state of rj in ^^O increases with 
the coupling constant. For 5^^^, > 0.90, however, the binding energy starts 
to decrease, and then ceases to exist when it becomes greater than 1.0. Fur- 
thermore, 3^(o^Ar) also decreases for 5^^^, > 0.85 and becomes repulsive at 
large values. 

The existence of at least one bound state of t] in medium mass nuclei of 
radius R = tqA^^^ can be understood by considering an equivalent square- 
well (complex) potential of depth Vq = Uq + iWo. One s— wave bound state 
is possible if the condition 

is satisfied [10] . In terms of the s-wave rjN scattering length a^^ , the con- 
dition is 

vr^ro A M, 



X<^ia^J<9X, X = ^[l + ^) , (5) 



and the potential is given by 



v,^-im.-4^)(i + ^)(^^L±J^]. (6) 
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Table 2. Parameters X and %X (both in fm) of the equivalent square- well potential 
calculated with = (0.28 + 0.19z) fm. 



Nucleus 


X 


9X 




0.109 


0.981 


160 


0.090 


0.810 


26 Mg 


0.065 


0.585 




0.048 


0.432 


902r 


0.029 


0.261 


208p|2 


0.016 


0.144 



In the above expression, the masses are in fm~ and Vq is in MeV. The values 
of X and 9X are given in table[2]for ro = 1.1 fm and a^^ = (0.28-|-0.19z) fm 
[3]. For nuclei with A < 10, one has to use the actual value of R which is 
substantially larger than the value given by I.IA^/^ fm. Additionally, the 
use of equivalent square-well potential and the omission of the imaginary 
part of the potential break down for lighter nuclei. Consequently, eqs.(l4])-([6]) 
should not be applied to nuclei with A < 10. Instead, detailed calculations 
should be performed. In Zr and Pb, 5R(a^^) > 9X. This explains why there 
are more than one bound state in these nuclei. 

2.2. Factorization Approximation 

In the factorization approximation (FA), the ryN scattering amplitude in 
eq. ([2]) is taken out of the Q integration and evaluated at a fixed momentum 
< Q > given by 

<'^>=- (^) m 

This choice of < Q > corresponds to a motionless target nucleon fixed 
before and after the ijN interaction. It has the virtue of preserving the 
symmetry of the t-matrix with respect to the interchange of k and k'. With 
this approximation, the ry-nucleus potential can be written as 

< k' I Vfa I k > = < k', -(k'+ < Q >) I t{\/^)nN^riN \ 

X |k,-(k+<Q>) >/(k'-k) , (8) 

where 

/(k' - k) = ^ I <p*{-k' - Q)0,(-k - Q) , (9) 
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Table 3. Binding energies and half- widths (both in MeV) of ?7-mesic nuclei obtained 
with the factorization approach for different values of the energy shift parameter 
A (in MeV). 



Nucleus 


Orbital (nf) 


A = 


A = 10 


A = 20 


A = 30 




Is 


-(2.18 + 9.96i) 


-(1.80 + 6.80i) 


-(1.42 + 5.19j) 


-(1.10 + 4.10i) 


16q 


Is 


-(4.61 + 11.57i) 


-(3.92 + 8. 13i) 


-(3.33 + 6. 37j) 


-(2.84 + 5. 17i) 




Is 


-(6.52 + 12.86i) 


-(5.63 + 9.16i) 


-(4.90 + 7.26j) 


-(4.29 + 5.96i) 




Is 


-(9.26 + 14.90i) 


-(8.09 + 10.75i) 


-(7.15 + 8.60j) 


-(6.37 + 7.13i) 




Is 


-(10.21 + 15. 4H) 


-(8.95 + 11.17i) 


-(7.94 + 8.97i) 


-(7.11 + 7.46i) 




Is 


-(10.84 + 15.70i) 


-(9.53 + 11.40i) 


-(8.49 + 9.18i) 


-(7.62 + 7.65i) 


32 g 


Is 


-(11.94 + 16.18i) 


-(10.56 + 11.80i) 


-(9.47 + 9.53j) 


-(8.55 + 7.97i) 


«Ca 


Is 


-(14.34 + 17.06i) 


-(12.75 + 12.55i) 


-(11.53 + 10.21i) 


-(10.51 + 8.59j) 




Is 


-(15.40 + 17.12j) 


-(13.73 + 12.66i) 


-(12.46 + 10.34i) 


-(11.40 + 8.73i) 




Is 


-(16.42 + 16.99i) 


-(14.65 + 12.63j) 


-(13.33 + 10.35i) 


-(12.24 + 8.77i) 




Ip 


-(2.10 + 15.09i) 


-(1.40 + 10.54i) 


-(0.75 + 8, 21i) 


-(0.19 + 6.62j) 




Is 


-(16.73 + 16.81i) 


-(14.94 + 12.53i) 


-(13.60 + 10.28i) 


-(12.51 + 8.72i) 




Ip 


-(2.68 + 15.00i) 


-(2.00 + 10.53i) 


-(1.34 + 8.25i) 


-(0.76 + 6.69j) 




Is 


-(17.04 + 16.88) 


-(15.23 + 12.59i) 


-(13.88 + 10.34i) 


-(12.77 + 8.78j) 




Ip 


-(3.17+ 14.82i) 


-(2.43 + 10.47i) 


-(1.73 + 8.24i) 


-(1.13 + 6.71i) 




Is 


-(18.57+ 17.34i) 


-(16.62 + 12.98j) 


-(15.20 + 10.69i) 


-(14.03 + 9.10i) 




Ip 


-(5.02 + 15.29i) 


-(4.20 + 10.94i) 


-(3.42 + 8.71i) 


-(2.75 + 7.18j) 




Is 


-(21.32 + 18.59i) 


-(19.15 + 13.97j) 


-(17.58 + 11.54i) 


-(16.29 + 9.84j) 


208 


Ip 


-(8.27+ 16.01i) 


-(7.19 + 11.47i) 


-(6.23 + 9.48i) 


-(5.40 + 7.94i) 


Is 


-(24.06 + 19.18i) 


-(21.88 + 14.44i) 


-(20.28 + 11.96) 


-(18.96 + 10.22i) 




2s 


-(4.89 + 11.04i) 


-(3.67 + 8.28i) 


-(2.81 + 6.79i) 


-(2.12 + 5.72i) 




Ip 


-(18.33 + 18.97i) 


-(16.31 + 14.27i) 


-(14.81 + 11.79i) 


-(13.56 + 10.06i) 




Id 


-(8.27+ 14.07i) 


-(6.17+ 10.56i) 


-(5.58 + 8.71i) 


-(4.66 + 7.41i) 



is the nuclear form factor having the normalization /(O) = A. Guided 
by the expression for the t-matrix in eq.([8]) is evaluated at = 

+ Mat — a = ^/sth — A, with A being an energy shift parameter. A 
downward shift of ~ 30 MeV is used to fit ttN scattering data [7j . 

The bound-state solutions obtained from using the covariant factorized 
potential with A = 0, 10, 20, 30 MeV are presented in table [3l The 
interaction parameters used in the FA calculations are same as those used for 
the off-shell calculations. The nuclear form factors used in the calculations 
can be found in refs.|lH[T^. A comparison between tables [T] and E] indicates 
that the FA results with A = 30 MeV are close to the off-shell results. This 
indicates that the ijN interaction in r] bound-state formation takes place at 
energies 30 MeV below the free-space threshold. 

The mass dependence of the Is binding energies for the A = 30 MeV 
case is shown in figured) The binding energies have been fitted empirically 
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1000 



Fig. 1 . Depenendence of the binding energy of ?7-mesic nucleus on the mass number 
A. 



with the formula 



aA + bA^/^ + cA^/^ + d, 12 < A < 208, 



(10) 



where a = 0.0003, b = -0.9562, c = 13.06, and d = -23.46 (all in MeV). 
The above equation resembles the semi-empirical mass formula with a as 
the volume energy term and b as the surface energy term. 



2.3. On-Shell Optical Potential 
First-order low energy r/-nucleus on-shell optical potential is 



< k' I Von I k >= 



1 



47r2/i 



/(k'-k). 



:ii) 



As will be seen later, Vqn corresponds to VpA with no energy shift (A = 0) 
and gives an upper limit to the value of Crj. The only input for the on-shell 
calculation is the s-wave rjN scattering length a^^. 

The scattering length is not directly measurable and its value is model 
dependent. Different models predict different values of o^^ [6]. The value 
of a ^ in the literature varies from 0.27 < 5R(a ^r) < 1.05, 0.19 < 9(a ^) < 
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Table 4. Binding energies and half-widths (both in MeV) of ry-mesic nuclei (Is 
state only) given by the on-shell optical potential for two different values of the 
scattering length a „ . No bound state exists in ''He. 



Nucleus 


F.A. (A = 0) 


a^j^ = (0.28 + 0.19i) fm 


a^^ = (0.51 + 0.21?) fm 


^He 


- 


- 


-(6.30 + 11.47f) 




— 


- 


-(3.47 + 6.79i) 


9Be 


— 


- 


-(13.78 + 12.45?) 






-(0.93 + 8.70) 


-(15.85 + 13.05?) 






-(2.71 + 10.9H) 


-(20.78 + 15.42?) 


12Q 


-(2.18 + 9.96i) 


-(2.91 + 10.22i) 


-(19.61 + 14.20?) 


160 


-(4.61 + 11.57i) 


-(5.42 + 11.43i) 


-(23.26 + 14.86i) 


20Ne 


-(6.52 + 12.86i) 


-(7.44 + 12.61i) 


-(26.72 + 15.94?) 




-(9.26 + 14.90i) 


-(10.34 + 14.40i) 


-(32.00 + 17.59?) 


26 Mg 


-(10.21 + 15.4H) 


-(11.24 + 14.76i) 


-(33.11 + 17.73?) 




-(10.84 + 15.70i) 


-(11.90 + 15.05i) 


-(34.06 + 17.96i) 


32s 


-(11.94 + 16.18i) 


-(13.10 + 15.6H) 


-(35.90 + 18.46?) 




-(14.34 + 17.06«) 


-(15.46 + 16.66?) 


-(38.85 + 19.16?) 




-(15.40 + 17.12?) 


-(16.50 + 17.07?) 


-(39.82 + 19.44?) 


52Cr 


-(16.42 + 16.99?) 


-(17.46 + 17.45?) 


-(40.5? + 19.76i) 


56 pg 


-(16.73 + 16.8H) 


-(17.76 + 17.52?) 


-(40.67 + 19.85?) 


58Ni 


-(17.04 + 16.88i) 


-(18.09 + 17.67i) 


-(41.13 + 20. 02i) 


70Zn 


-(18.57 + 17.34^;) 


-(19.63 + 18.47i) 


-(43.52 + 20.93i) 


90 2p 


-(21.32 + 18.59i) 


-(22.41 + 19.97i) 


-(48.40 + 22. 60i) 


208pb 


-(24.06 + 19.18i) 


-(24.55 + 19.57i) 


-(50.27 + 21.42?) 



0.37. A reason for this large range is unavailability of data on rjN elastic 
scattering, which is essential in determining the value of a^^ . 

The binding energies and half- widths given by eq. ()lip are presented in 
table m for two different values of a^^ . The nuclear form factors are the 
same as those used in the FA calculations. For these two scattering lengths, 
no bound state can exist in ^He. Upon comparing the third column of 
table m with the off-shell calculation (table [T|) , it can be seen that the on- 
shell approximation predicts more strongly bound ?7-mesic nuclei. Also, as 
expected, the on-shell results for a^^ = (0.28 + 0.19?) fm are similar to those 
of FA with A = MeV. 
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Fig. 2. The calculated energy dependences of the DCX reaction described in the 
text leading to the double isobaric analog state. 

3. Eta-Mesic Nucleus and Pion DCX 

Eta-mesic nucleus can affect high-energy pion double-charge-exchange 
(DCX) reactions through the formation of "mesonic compound nucleus" [13j . 
For pion kinetic energies > 400 MeV, r] production channel is open in 
most nuclei. As a result, the DCX reaction can proceed via vr"'" — > vr'' — > vr" 
or vr"^ — > r/ ^ 7r~. While 7r° is in the continuum, r/ can either be in the con- 
tinuum or in a nuclear bound state. The DCX amplitudes associated with 
the ry-nucleus bound states, on the other hand, have resonance structure. 

The calculated cross sections for the reaction ^^C(7r^, 7r~)^^N as a func- 
tion of Tjr at momentum transfers g = and 210 MeV/c are shown in figure 
[2j The solid curves represent the contribution from the resonant ampli- 
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tude associated with the formation of ry-mesic nucleus. The dashed curves 
represent contribution from the nonresonant amphtude only. The interfer- 
ence of these amplitudes is responsible for the presence of narrow resonance 
structure at ~ 415 MeV for q = 210 MeV/c case. The width of the 
resonance is about 10 MeV, which reflects the width (~ 11 MeV) of the 
77-mesic nucleus used in the calculations. The DCX studies, therefore, can 
be used as an alternative way to determine the width of r/-mesic nucleus. 
Other processes in which one can expect to see these kind of effects is (vr, tt') 
reactions leading to certain specific final states. The r/-nucleus bound state 
amplitude for these kind of reactions is not small, in comparison to the 
nonresonant amplitude. The study of the resonance pattern in the energy 
dependence of the cross section can yield information on the relative phase 
between resonant and nonresonant amplitudes. 

4. Experimental Search for Eta-Mesic Nucleus 

Several experiments were performed to detect 77-mesic nucleus. The 
first experiment at Brookhaven National Laboratory [Hj in 1987 was based 
on the work of Liu and Haider [T5|. A second experiment, based on the 
pion DCX calculations [13] was done at LAMPF by Johnson et al. |16| . 
While these experiments could not confirm unambiguously the existence of 
ry-mesic nucleus, they did not rule out such a possibility either. In recent 
publications, Sokol et al. [T7] claim to have observed ?7-mesic nucleus in 
experiments involving photo-mesonic reactions. 

Recently, an experiment has been done at Jiilich by the COSY-GEM 
Collaboration [18], making use of the transfer reaction p-|-^^Al— >^He-|-^^Mg^. 
Analysis of the data clearly indicates the detection of bound state of r/ in 
^^Mg. Vigorous efforts are also underway to detect bound state of 77 in 
^'^He. It is claimed that a "quasibound state" of rj in ^He may have been 
observed in the dp ^^Her/ reaction [191 120] . This is in sharp contrast to the 
findings of the COSY-11 and COSY-at-WASA groups [H], also looking for 
the bound state of rj in ^He in dp collisions. Within the statistical sensitiv- 
ity achieved by them, they could not confirm the existence of 77— ^He bound 
state. This is in agreement with the predictions presented in this paper. 

5. Conclusions 

The formation of r/-mesic nucleus is a natural consequence of the attrac- 
tion between the i] meson and the nucleon at very low energies. However, 
the attraction is not strong enough to bind an 77 onto a single nucleon. The 
bound state formation is possible in a finite nuclei with mass number greater 
than 10. 
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The calculated binding energies and widths of r/-nucleus bound states 
strongly depend on the subthreshold dynamics of the rjN interaction. The 
present analysis indicates that the average r]N interaction energy in mesic- 
nucleus formation is below the threshold. What matters for the bound-state 
formation is not the rjN interaction at the threshold but the effective in- 
medium interaction. Because the subthreshold behavior of rjN interaction 
is very model dependent, it is useful for theorists to publish not only the 
r/-nucleon scattering length a^^^ but also the corresponding subthreshold 
values as a function of the shift parameter A. 

The downward shift in the effective interaction energy can lead to a 
substantial reduction of the attraction of in-medium r/-nucleon interaction 
with respect to its free-space value. Consequently, predictions based upon 
using free-space rjN scattering length inevitably overestimate the binding 
of rj. One must bear this in mind when using the predictions given by such 
calculations as guide in searching for ry-nucleus bound states. 

Recent experimental confirmation |18| of the existence of r/-mesic nucleus 
will enable us to improve upon the existing models or identify additional 
physics that has to be incorporated in them. Because the binding energies of 
r] meson depend strongly on the coupling between the r]N and the A^*(1535) 
channels [3] , studies of ?7-mesic nucleus will be able to yield detailed informa- 
tion on the rjNN* coupling constant involving bound nucleons. It can also 
lead to a new class of nuclear phenomenon, ?7-mesic compound nucleus res- 
onances. An awareness of this phenomenon could be helpful to the analysis 
of nuclear reactions at energies above the threshold for r] production. 

The ?7-mesic nuclear levels correspond to an excitation energy of ~ 540 
MeV, to be compared with an excitation energy of ~ 200 MeV associated 
with the A- and S-hypernuclei. The existence of nuclear bound states with 
such high excitation energies provides the possibility of studying nuclear 
structure far from equilibrium. 
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